Modified spin-wave description of the nuclear spin relaxation in ferrimagnetic 

Heisenberg chains 
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We make a modified spin-wave description of the nuclear spin relaxation in Heisenberg 
alternating-spin chains with antiferromagnetic exchange coupling. In contrast with the conven- 
tional one- dimensional antiferromagnetic spin-wave theory, which is plagued with the divergence of 
the sublattice magnetization even in the ground state, the present spin-wave description is highly 
successful over a wide temperature range. The temperature dependence of the relaxation rate Tj" 
significantly varies with the crystalline structure, exhibiting both ferromagnetic and antiferromag- 
netic aspects. T^ 1 further shows a unique dependence on the applied field, which turns out an 
indirect observation of the quadratic dispersion relations. 

PACS numbers: 75.10.Jm, 75.40. Mg, 75.30. Ds, 76.60. Es 



I. INTRODUCTION 

Quantum mixed-spin chains with magnetic ground 
states, namely, quantum ferrimagnets, are one of the hot 
topics and recent progress |l|-|25t in the theoretical under- 
standing of them deserves special mention. Coexistent 
ferromagnetic and antiferromagnetic long-range orders 
P] in the ferrimagnetic ground state in particular interest 
us. The ground-state magnetizations of antifcrromagnets 
and ferromagnets are zero and saturated, respectively, 
and therefore ferrimagnets may be recognized to possess 
in-between ground states. Hence the ground-state excita- 
tions in ferrimagnets are twofold JlC],[llJ • The elementary 
excitations of ferromagnetic aspect, reducing the ground- 
state magnetization, form a gapless dispersion relation, 
while those of antiferromagnetic aspect, enhancing the 
ground-state magnetization, are gapped from the ground 
state. The dual structure of the excitations results in 
unique thermal behaviors : The spe- 

cific heat and the magnetic susceptibility times temper- 
ature behave like T 1 ' 2 and T _1 at low temperatures, re- 
spectively, whereas they exhibit a Schottky-like peak and 
a round minimum at intermediate temperatures, respec- 
tively. Quantum ferrimagnets in a magnetic field provide 
further interesting issues such as the double-peak struc- 
ture of the specific heat and quantized plateaux 
in the ground-state magnetization curves [T^JTf^ , |20| , ^l[| 
In particular it has quite recently been reported pa ] 
that ferrimagnetic Heisenberg chains can exhibit multi- 
plateau magnetization curves even at the most symmet- 
ric point, that is, without any anisotropy and any bond 
polymerization. 

It is true that theoretical investigations into quantum 
ferrimagnets are now active and interesting in them- 
selves, but we should still be reminded that such vig- 
orous arguments more or less originate in the pioneer- 
ing efforts |26|j2^ ] to synthesize bimetallic materials in- 
cluding one-dimensional systems. The first ferrimag- 
netic chain compound @, MnCu(dto)2(H 2 0)3-4.5H 2 
(dto = dithiooxalato = S2C2O2), was synthesized by 



Gleizes and Verdaguer and stimulated the public interest 
in this potential subject. The following examples [^9|j30| 
of an ordered bimetallic chain, MnCu(pba)(H 2 0)3-2H20 
(pba = 1, 3-propylenebis(oxamato) = C7H 6 N 2 06) 
and MnCu(pbaOH)(H 2 0) 3 (pbaOH = 2-hydroxy-l, 3- 
propylenebis(oxamato) = C7H6N2O7), exhibiting more 
pronounced one dimensionality, activated further physi- 
cal [pT|-^3|, as well as chemical |p4f , investigations. The 
serial chemical explorations condensed into the crystal 
engineering of a molecule-based ferromagnet |35|— the as- 
sembly of the highly magnetic molecular entities within 
the crystal lattice in a ferromagnetic fashion. 

Thus, a good amount of chemical knowledge on quasi- 
one-dimensional quantum ferrimagnets has been accumu- 
lated and static properties of them have been revealed 
well. However, little is known about dynamic proper- 
ties of quantum ferrimagnets. To the best of our knowl- 
edge, in the theoretical field, it was not until quite re- 
cently that the dynamic structure factors were calculated 
jnj, while in the experimental field, any direct observa- 
tion of the energy structure is not yet so successful, for 
instance, as that J36| for the Haldane antiferromagnets 
p7| , j38|| . In such circumstances, Fujiwara and Hagiwara 
performed |39| nuclear-magnetic-resonance (NMR) mea- 
surements on bimetallic chain compounds. The measured 
temperature and applied-field ranges were rather limited 
and their argument was not so conclusive. However, they 
suggested that the nuclear spin relaxation could be a 
useful tool in order to look into the low-energy struc- 
ture peculiar to quantum ferrimagnets. In response to 
this stimulative experiment, here we calculate the nuclear 
spin relaxation rate T^ 1 in terms of a modified spin- wave 
theory and strongly encourage further experimental in- 
vestigations. 



II. FORMULATION 



We describe alternating-spin chain compounds by the 
Hamiltonian 
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3=1 i=i 

(2.1) 

where S- = 5(5 + 1), s| = s(s + 1), and we have 
set their g factors both equal to g because the dif- 
ference between them amounts to at most several per 
cent of themselves in practice p()| . We further set 
the unit-cell length, which is twice the lattice constant, 
equal to unity in the following for the convenience of 
calculation. Magnetic properties of the ferrimagnetic 
family compounds such as MCu(pba)(H 2 0)3-2H 2 and 
MCu(pbaOH)(H 2 0) 3 (M = Mn, Ni) are well described 
within this isotropic Hamiltonian p9pc| , |40| , fll| . Consid- 
ering the electronic-nuclear energy-conservation require- 
ment, the direct (single-magnon) process is of little sig- 
nificance but the Raman (two-magnon) process plays a 
leading role in the nuclear spin-lattice relaxation p2[ . 
Then the relaxation rate is generally given by 



J_ 47r(g/i B ^7N) 2 



-E n /k B T 



x \{m\Y,,(A z j S] + a z j s z j )\n)\'S(E m -E n -nLJ U ), (2.2) 

where A* and o| are the dipolar coupling constants be- 
tween the nuclear and electronic spins in the jth unit cell, 
cjn = InH is the Larmor frequency of the nuclei with 7^ 
being the gyromagnetic ratio, and the summation J2 n ^ s 
taken over all the electronic eigenstates \n) with energy 
E n . 



In order to rewrite the Hamiltonian (2.1) within the 
framework of the spin wave theory, we introduce the 
bosonic operators for the spin deviation in each sublattice 
via 
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b}j2s-b} bj , 



Sj — - S djQjj , 

8»j =-s+ b]b 3 , 



(2.3) 



where we regard 5 and s as quantities of the same order. 
Now we obtain the bosonic Hamiltonian as 



HsW = -Eclass + Ho + Hi + 0(S 1 ) , 



(2.4) 



where B c iass = —2sSJN is the classical ground-state en- 
ergy, and Ho and Hi are the one-body and two-body 
terms of the order 0(S 1 ) and 0(S ), respectively. We 
may consider the simultaneous diagonalization of Ho and 
Hi in the naivest attempt to go beyond the linear spin- 
wave theory. However, such an idea ends in failure bring- 
ing a gap to the lowest-lying ferromagnetic excitation 
branch. Thus we take an alternative approach |l3| at 
the idea of first diagonalizing Ho and next extracting 
relevant corrections from Hi . Ho is diagonalized as 



Ho 



En 



E 



(2.5) 



where E — J^ki^k — + s )l i s tne ^(S 1 ) quantum 
correction to the ground-state energy, and a\ and f3 k are 
the creation operators of the ferromagnetic and antifer- 
romagnetic spin waves of momentum k whose dispersion 
relations are given by 



uj k ± (5 - s)J =f gm H , 



with 



oJk = JJ(S - s) 2 +4Sssin z (k/2). 



(2.6) 



(2.7) 




FIG. 1. Dispersion relations of the ferromagnetic and 
antiferromagnetic elementary excitations, namely, the low- 
est-energy states in the subspaces of M = N/2 =p 1. The 
linear- and interacting-spin-wave calculations are shown by 
the dotted and solid lines, respectively, whereas x represents 
the quantum Monte Carlo estimates (N — 32). Here we plot 
the excitation energy E(k) taking the ground-state energy as 
zero. 
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The Wick theorem allows us to rewrite Hi as 



Hi = Eh. - £ (5^a{a k + duj+P^k 

k 



(2.8) 



where iArei contains irrelevant terms such as ctkflk and 
7i quar t contains residual two-body interactions, both of 
which are neglected in the following so as to keep the 
low-energy structure qualitatively unchanged. E\ — 



-2jN[r{ 



^/S)r 1 r 2 ] is the O(S ) cor- 



rection to the ground-state energy, while 



2(5 + s)r 3 



sin 2 (fc/2) 



A 



K ±(s- s)} 



(2.9) 



are those to the dispersions, where the key constants A 
and A are defined as A = (27V)- 1 ^[(S + s)/w fc - f] 
and A = -iV" 1 X)fe(v / 5s/w* : )cos 2 (A;/2). The resultant 
Hamiltonian is compactly represented as 
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sw 
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(2.10) 
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with Eg = Aiass + E a + Ei and cD fc = u _ 

We show in Fig. [l] the linear- and interacting-spin-wave 
dispersions, and aJjr, together with the numerical so- 
lutions |^| obtained through imaginary-time quantum 
Monte Carlo calculations p3-45|. The spin- wave descrip- 
tion of the low-energy structure is fairly good. Even the 
linear spin waves allow us to have a qualitative view of 
the elementary excitations. The relatively poor descrip- 
tion of the antiferromagnetic branch by the linear spin 
waves reminds us of the spin-wave treatment of mono- 
spin Heisenberg chains, where the theory accurately de- 
scribes ferromagnetic chains, while it only gives a qual- 
itative view of antiferromagnetic chains. The spin-wave 
approach to the present system is highly successful any- 
way for both excitation branches. The spin-wave series 
potentially lead to the goal even for the antiferromagnetic 
branch. The high applicability essentially originates in 
the fact that the spin deviations 
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dfc, (2.11) 



with ( ) g denoting the ground-state average, no more 
diverge in the present system with S ^ s. We are con- 
vinced that the quantity A should be recognized as the 
quantum spin reduction. 

In terms of the spin waves, the relaxation rate ( [2.2] ) is 
expressed as 



t, = -^?(3mb7n) 2 E E s ^+ q 

k,q o=± 

x [(^cosh6> fc+9 cosh6> fc ) 2 <« +(? + 1) 
+(a z q sivhe k+q sivhe k ) 2 n% +q (nl + 1) 
-2A z q a z q (cosh9 k smhe k ) 2 nl(nl + 1)] , 



(2.12) 



where n k = (u\a.k) and n k = (P^Pk) are the thermal av- 
erages of the boson numbers at a given temperature, and 

A q = J2j eiqU ~ 1/4)A j aild a q = J2j eiqU + 1/4)a j arC tllC 

Fourier components of the hyperfine coupling constants. 
Taking into account the significant difference between the 
electr onic and nuclear energy scales (hui^ < 10~ 5 J), Eq. 
(|2.12|) ends in 
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x [(A 2 cosh 2 6» fc - a*sinh 2 fe ) 2 n fe (n k + 1) 
+ (A z smh 2 9 k - a z cosh 2 9 k ) 2 n+(n+ + 1)] , (2.13) 



where we have assumed little fc-dependence of A z q and a z , 
and thus replaced A; l=_ 2 k ano - a q=-2k A q=o = A * and 
a q=Q = a z , respectively. 

The estimation of the relaxation rate is now reduced 
to the calculation of the spin-wave distribution functions. 
Though the ground-state distribution is well controlled, 
the naivest thermodynamics, based on the partition func- 
tion Z = Tr[e~ Wsw / fcBT ], breaks down as temperature 
increases. Hence we modify the spin-wave theory j49[ ] 
introducing an additional constraint in minimizing the 
free energy. Requiring the total magnetization to be 
zero, Takahashi |46| obtained an excellent description of 
the low-tempcraturc thermodynamics of one-dimensional 
Heisenberg ferromagnets. Taking his core idea but re- 
placing the ferromagnetic constraint + s j) = 
by 



E 

3 



({S z - s z ) + 2A) = (S + s) (N - £ £ £ = 



(2.14) 



we obtain the modified spin-wave distribution functions 
as 



1 



(2.15) 



with a Lagrange multiplier fi due to the condition (2.14). 
In comparison with the ferromagnetic cases, the quantum 
correction 2A is necessary for ferrimagnets. The spin- 



■ s under ferromag- 

■ = S+s — 2 A un- 



wave treatment gives {Si- + s| ) g = S 
netic exchange coupling but \S Z — s z 
der antiferromagnetic exchange coupling. Indeed, with- 
out the quantum correction, we reach a quite poor de- 
scription of the thermal quantities O] . We show in Fig. 
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g the modified spin-wave calculations of the magnetic- 
susceptibility-temperature product, which is closely re- 
lated with the relaxation rate ]5Q] and is given in terms 



of nf as 



XT 



EE 

k <J=± 



<« + !)■ 



(2.16) 



The obtained results are fairly successful considering 
that these are the spin-wave calculations in one dimen- 
sion. The modified spin waves well reproduce the low- 
tcmpcraturc ferromagnetic divergence, which is propor- 
tional to T , the high-temperature antiferromagnetic 
increase toward the paramagnetic behavior [S(S + 1) + 
,s(s+l)]/3, and therefore, the round minimum at interme- 
diate temperatures. In particular, the low-temperature 
description by the interacting spin waves may be re- 
garded as accurate. Now we proceed to the argument of 
the relaxation rate in terms of the modified interacting- 
spin-wave theory. Since the applied field H is so small in 
practice ^] as to satisfy g/iBH <10_ 2 J, we neglect the 
Zeeman term of the Hamiltonian ( [2.1C| ) in the estimation 
of nf. 
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FIG. 2. Temperature dependences of the zero-field mag- 
netic susceptibility times temperature in the case of 
(S, s) — (1, |). x represents the quantum Monte Carlo es- 
timates, whereas the dotted and solid lines show the modi- 
fied-spin-wave calculations starting from the linear- and in- 
teracting-spin-wave dispersion relations, respectively. The 
low-temperature behaviors are scaled up in the inset. 
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FIG. 3. Temperature dependences of the nuclear spin re- 
laxation rate at various values of the applied magnetic field 
in the case of (S, s) = (1, |). 
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FIG. 4. The momentum distribution functions of the 
ferromagnetic (a) and antiferromagnetic (b) spin waves 
in the modified spin-wave theory starting from the in- 
teracting-spin-wave dispersion relations in the case of 
(5, S ) = (1,|). 



III. RESULTS 

Here still remains an adjustable parameter A z /a z = r. 
The dipolar coupling is quite sensitive to the location of 
the nuclei because the coupling strength is proportional 
to g?~ 3 with d being the distance between the interacting 
nuclear and electronic spins. In the proton-NMR mea- 
surements on NiCu(pba)(H20)3-2H20 for instance, 
it was demonstrated that the protons relevant to the re- 
laxation rate do not originate in the H2O molecules but 



lie in the pba groups which are located beside the Cu ions. 
Thus, for these family compounds, r may reasonably be 
set equal to zero. We show in Fig. || the corresponding 
calculations. As the measurements were performed at 
rather high temperatures, which are beyond the quantita- 
tive reliability of the present theory, it is impossible to fit 
the calculations to the experimental findings. However, 
the calculations at r = well explain the observations of 
monotonically decreasing behaviors of T^ 1 as a function 
of temperature, which are in contra st wi th the app ear- 
ance of %T. When we compare Eqs. ( 2.13| ) and ( |2.16| ), we 
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realize that T^ 1 could tell more than \T due to its ad- 
justable prefactors to n£(7i£ + 1). Setting r equal to 0.4, 
we obtain temperature dependences of T^ 1 exhibiting a 
round minimum, where the antiferromagnetic spin-wave 
contribution is much more accentuated than the ferro- 
magnetic one nZ. Let us observe the momentum depen- 
dences of in Fig. ^. nZ exhibits a sharp peak around 
k = at low temperatures, which is rapidly reduced with 



600 



+ 



is 



the increase of temperature. On the other hand, n 
an increasing function of temperature at an arbitrary mo- 
mentum, though its broad peak is smeared out with the 
increase of te mperature. The field-dependent prefactor 
in Eq. ( 2.1 3| ) , coming from the ene rgy conservation re- 
quirement S(E m — E n — hutti) in Eq. ( |2.2|) , predominantly 
extracts the k = components from n k . Therefore, the 
parametrization r — 0.4 ~ | = tanh 2 #fc = o results in the 
strong suppression of the ferromagnetic aspect of T-f . 
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FIG. 5. Field dependences of the nuclear spin relax- 
ation rate at various values of temperature in the case of 
= (l,i). 
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FIG. 6. Temperature dependences of the nuclear spin relaxation rate at various values of the applied magnetic field in various 
constituent-spin cases. 



Another motivation of the NMR measurements on fer- 
rimagnetic chains may be the characteristic field depen- 
dences of T x _1 shown in Fig. ||. It is due to the quadratic 



dispersion relations (2.6) that T x depends on the ap- 



plied field. Hence the present field dependence allows 
us to look into the low-energy structure peculiar to fer- 
rimagnets. The high linearity of T^ 1 with respect to 
H^ 1 / 2 denotes the predominance of the k ~ compo- 
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nents in the fc-summation in Eq. (2.13). The predomi- 
nance is reduced with the increase of H and finally there 
arises a logarithmic field dependence of T^ 1 from the k- 
integration of [(Ssk) 2 + 2(S-s)Sshuj N /J}- 1 ^. If we take 
J/k B = 121[K] @ relevant to NiCu(pba)(H 2 0) 3 -2H 2 0, 
Fig. H suggests that the logarithmic behavior should ap- 
pear under H > 10[T]. The T\ measurements ]3!| on 
NiCu(pba)(H 2 0) 3 -2H 2 at 280[K] with an applied field 
up to 3.15[T] resulted in a monotonic linear dependence 
of Tj on H^ 1 ! 2 . The high-field measurements at lower 
temperatures are expected. 

In order to stimulate extensive NMR measurements on 
ferrimagnetic compounds, we show in Fig. ^temperature 
dependences of T{~ in general constituent-spin cases. We 
select three particular values of r: (a) r = 0, where the 
relevant nuclei are located closely to the smaller spins s 
and the ferromagnetic and antiferromagnetic mixed na- 
ture is displayed; (b) r = tanh 0fe = o, where the rele- 
vant nuclei are located near the smaller spins s rather 
than the larger spins S and the ferromagnetic aspect 
is strongly suppressed; and (c) r = coth 0fe=o, where 
the relevant nuclei are located near the larger spins S 
rather than the smaller spins s and the antiferromag- 
netic aspect is strongly suppressed. Thus, the chemi- 
cal technique might enable us to extract separately the 
ferromagnetic/antiferromagnetic feature of ferrimagnets 
from Tj~ . Let us observe the temperature dependences 
from the point of view of the low-energy structure, which 
has been revealed in Fig. |l|. The distinct behaviors 
Fig. ^(b) and Fig. ^|(c) are reminiscent of the antifer- 
romagnetic and ferromagnetic yT products, respectively. 
The antiferromagnetic aspect should be dominated by 
the antiferromagnetic gap uj~£_ . Although the spin- wave 
description within the up-to-0(S' 1 ) approximation con- 
siderably underestimates the antiferromagnetic gap, its 
estimate wt_ = 2(5 — s) J can be a qualitative guide. 
We learn that the antiferromagnetic gap is in proportion 
to S — s and indeed find the slower activation and the 
higher-located peak of the antiferromagnetic component 
of T^ 1 for (S, s) — (|, |) in comparison with those for 
(S, s) — (1, |) and (5, s) = (|, 1). A careful observation 
of Fig. |6](a) shows that T^ 1 reaches a minimum around 
fesT/J = 2.5 for (5, s) — (|, In the two other cases, 
T^ 1 is still on the decrease at k B T/J = 6.0. We may 
understand whether T^ 1 exhibits a minimum— that is to 
say, a ferromagnetic-to-antiferromagnetic crossover— at a 
tractable temperature in connection with the balance be- 
tween the ferromagnetic band width W~ and the an- 
tiferromagnetic gap A. The ferromagnetic decreasing 
tail exists for k B T < W~ , whereas the antiferromag- 
netic increasing behavior is remarkable for k B T < A. 
If we evaluate W~ and A as ui 



k—TT 



UJ 



2Js and 

wjjL-o = 2(5 — s) J, respectively, we may expect a de- 
tectable minimum of T^ 1 as a function of temperature 
for 2s < 5. The larger S — 2s, the more pronounced 
crossover may be detected. 



IV. SUMMARY AND DISCUSSION 



A relaxation mechanism based on the interaction with 
spin waves has so far been applied to magnetic insula- 
tors only in a temperature range far below the onset 
of the long-range order p2[ . Wc again stress that the 
present argument is a positive use of the spin-wave the- 
ory in one dimension in a temperature range above the 
onset of the (three-dimensional) long-range order. The 
quantum divergence of the spin reduction, inherent in 
one-dimensional antiferromagnets, no more plagues the 
present system, while we have avoided the thermal di- 
vergence of the bosonic distribution function modifying 
the spin-wave theory, that is, imposing a certain con- 
straint on the magnetization. Though the present anal- 
ysis must be a qualitative guide for the experiments over 
a wide temperature range, the description may be very 
precise especially for k B T/ J < 0.2. T\ measurements at 
low temperatures and/or under high fields are strongly 
encouraged. 

The experimental development depends on the syn- 
thesis of relevant materials. Although it is the pioneer- 
ing measurements |m| on NiCu(pba)(H 2 0)3-2H 2 that 
have motivated the present study, the family compounds 
MCu(pba)(H 2 0)3-2H 2 may not be so useful as to ver- 
ify the present analysis. The linewidth in the NMR spec- 
tra for NiCu(pba)(H 2 0)3-2H 2 considerably broadens at 
low temperatures and therefore the extraction of T\ was 
restricted to a temperature range k B T/J > 0.5, where 
it is rather hard for the present theory to make a quan- 
titative interpretation, namely, a precise determination 
of the geometric parameters A z and a z . Even though 
the linewidth broadening is inevitably dominated by the 
crystalline structure, a strong exchange coupling must be 
desirable here anyway. The idea of designing ligands ca- 
pable of binding to two different metal ions with different 
donor atoms has not yet ended in an exch ang e coupling 
constant J beyond 85[cm _1 ] ~ 122fc B [K] [f27||. A break- 
through may be made by bringing into interaction metal 
ions and stable organic radicals. Ceneschi et al. (si] |5j| 
synthesized a series of ferrimagnetic chain compounds of 
general formula M(hfac) 2 NITR, where metal ion com- 
plexes M(hfac) 2 with hfac = hcxafluoroacetylacetonate 
are bridged by nitronyl nitroxide radicals NITR. Their 
exchange coupling constants significantly vary with M 
and R but are in general fairly large in comparison with 
those of bimetallic chain compounds. An antiferromag- 
netic coupling of 313[cm _1 ] ~ 450fcB[K] was obtained for 
M = Mn (S = |) and R = isopropyl (s = i) @], while 
that of 4241cm- 1 ] ~ 610fc B [K] for M = Ni (5 = 1) and 
R = methyl (s = i) |54| . Neither compound shows any 
transition to the three-dimensional magnetic order down 
to k B T I J = 0.02, which is also suitable for our purpose, 
apart from designing a molecule-based ferromagnet. We 
theoreticians hope as many experimentalists as possible 
will take part in this exciting business— dynamic proper- 
ties of one-dimensional quantum ferrimagnets. 
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